Quantum Spinon Oscillations 
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The full quantum dynamics of a spinon under external magnetic fields is investigated by using the 
time-evolving block decimation (TEBD) method within the microcanonical picture of transport. 
We show that the center of the spinon oscillates back and forth in the absence of dissipation. 
The quantum many-body behavior can be understood in a single-particle picture of transport and 
Bloch oscillations, where quantum fluctuations induce finite life times. Transport, oscillations and 
lifetimes can be tuned to some degree separately by external fields. Other nontrivial dynamics such 
as resonance as well as chaos have also been discussed. 
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Controlling the time-evolution of quantum states and 
manipulating the dynamics of quantum matter have at- 
tracted considerable theoretical and experimental atten- 
tion in recent years due to their relevance to fundamental 
physics as well as potential applications in information 
storage, encoding and processing. Up to now, most ef- 
forts in this field are based on solid state systems. Re- 
cently, the progress in ultracold atomic gases in opti- 
cal lattices not only provides an ideal platform for simu- 
lating quantum many-body models in condensed matter 
physics, but also paves the way for quantum manipu- 
lation. Due to the unique properties of the ultracold 
atomic gases in optical lattices, such as the extremely 
low dissipation rates and long coherence times, we can 
explore the dynamics of quantum many-body systems in 
a clean and dissipation- less environment. This enables 
us to manipulate the quantum many-body state with an 
unprecedented degree of precision, and study new phys- 
ical effects [ij. An early breakthrough in the field was 
the first direct observation of Bloch oscillations (BO) in 
a tilted optical lattice 0, In solid state systems, scat- 
tering by impurities would result in strong damping of 
Bloch oscillations, while in an optical lattice the perfect 
optical crystals are free from any imperfections, and what 
is more, the long coherence times and high tunability of 
optical lattices not only enable us to observe BO directly, 
but also to manipulate the dynamics of the quantum par- 
ticles via external fields. 

Compared with the dynamics of the quantum particles, 
it may be more interesting to explore and manipulate the 
quasi-particles, which emerge as elementary excitations 
in quantum systems consisting of a collection of interact- 
ing particles, and may behave differently compared to the 
original particles. The most interesting example of this 
is known as fractionalization: the particles are effectively 
split into smaller constituent quasi-particles, which only 
carry a fraction of quantum numbers. 



In this paper, we study the quantum dynamics of one 
of the most known quasi-particles: the spinon, which is 
an excitation separating two degenerate states of oppo- 
site magnetization in a quantum spin chain and usually 
carries a fractionalized quantum number (spin-^). Using 
the time-evolving block decimation (TEBD) methodN, |5| 
within the microcanonical picture of transport ^6^, '7f7we 
study the time evolution from a quantum many-body ini- 
tial state with a spinon, and show how to manipulate its 
dynamics via external fields. In particular, we find that 
this quantum many-body problem can be interpreted in 
the picture of a quantum quasiparticle at the single- 
particle level, which shows controllable linear transport 
and (Bloch) oscillation behavior depending on perpen- 
dicular external fields, where the field along the magne- 
tization direction controls the Bloch oscillation physics 
and the transverse field controls the linear velocity. The 
broadening of the quasiparticle is also tunable by the field 
strength. We also find a quantum resonance behavior 
of the spinon under a periodically driven external field. 
This may pave the way to the controlled manipulation 
of quasiparticles, also in view of potential applications in 
condensed matter and ultracold atom physics. 

Our Hamiltonian is a finite ID transverse Ising model 
with an additional magnetic field along the z-direction: 
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(1) 



with a ferromagnetic coupling J — —1, Pauli matri- 
ces cr^'^ on sites i, and a transverse field along the x- 
direction that breaks magnetization conservation. With- 
out the magnetic field haf, the transverse Ising model 
can be solved exactly by a Jordan- Wigner transforma- 
tion and is known as a classic paradigm of a quantum 
phase transition[8]. 

Recently, the transverse Ising model has been realized 
experimentally in cobalt niobate and a Eg, symmetry has 
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FIG. 1: The time evolution of the spinon for h — and (a) 
3 = 0.1, (b) 5 = 0.3. 



been observed in the vicinity of the quantum critical 
point of this model In the presence of a mag- 

netic field along the z-direction h, Eq.([T]) can no longer 
be solved exactly by Jordan- Wigner transformation, be- 
cause it would induce a nonlocal term. In this paper, h 
is a constant or a time-dependent field, and in both case, 
it will lead to nontrivial dynamics of the spinon. 
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FIG. 2: (a) Magnetization Sf at different times (black 
squares: t = 0.05; red dots: t = 10; green upward trian- 
gles: t = 20; blue downward triangles: t = 30) for h — 0, 
g — 0.3. (b) The dependence of the spinon velocity on g. 



Our many-body initial state is prepared as follows: 
a spinon with a finite velocity is located at the center 
of lattice. To realize this initial state, we can choose 
its wavefunction as \ip)t=o = 7^ Aje''^°*|z), where A^o 
is the normalization constant and \i) denotes a perfect 
static spinon located between site i and i + 1: |j) = | t 
)i • • • I t)i| • • • I 4-)l, where L is the length of the 

lattice. If we choose Ai — 1, the initial state is a per- 
fect Bloch wave function with definite wave vector fco. 
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FIG. 3: The time evolution of a spinon for g = 0.3 and (a) 
h = 0.1, (b) h = 0.2. 



and the spinon is totally delocalized over the length of 
the lattice. To ensure that the spinon is initially located 
at the center of the lattice, we then choose a Gaussian 



A,. 



-2(i-L/2Y 



Given the exponentially fast decay of 



Ai away from the center of the lattice, we can just keep 
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L/2 



Al/2+1, Al/2- 



and make all other Ai = for con- 



structing the initial wave function. In our calculation we 
furthermore choose ko — Tr/2. We also define the average 
position of the spinon P as the position where Sz = 0. 
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FIG. 4: The time evolution of the position of the spinon for 
ft = 0.2 and different g. 

As certain excitations can be modeled as quasiparti- 
cles, i.e., having similar properties to real particles, we 
may conjecture that this is the case here and compare 
the spinon to a real particle, where the field coupling to 
lead to the motion of the particle and where the field 
coupling to az acts as a constant electric field along the 
axis in the particle picture, as the Zeeman energy induced 
by haz is proportional to the displacement of the spinon 
from its original position. To verify this picture quanti- 
tatively, we use TEBD within the microcanonical picture 
of transport to calculate the time evolution of the spinon 
from our initial state \(p)t=o which is not an eigenstate 
of the Hamiltonian. In the course of the real time evolu- 
tion we use the length of our lattice L = 30 and take the 
truncation dimension ^ = 80 and time step At = 0.05. 
The convergence is checked by taking larger x- Times are 
measured in inverse hopping strengths \J\~^ throughout 
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the paper. 




FIG. 5: (a) The time evolution of the position of tlie spinon 
for g = 0.3 and different h; (b) the hnear relation between the 
frequency of BO and h. 

First we focus on the case h = 0; the real-time evolu- 
tion from our initial state \ip)t=o for different parameters 
g is shown in Fig. [TJ (a) g = 0.1; (h) g — 0.2. We find 
that the spinon propagates along the — x-direction with 
a velocity proportional to g as shown in Fig. [5] (a) . The 
linear relation between the velocity and g has been nu- 
merically verified in Fig. [2{b) . The quantum fluctuations 
increase the width of the spinon during the time evolu- 
tion, giving a lifetime to the quasiparticle. The quantum 
fluctuations make the quantum spinon different from its 
classical counterpart: the magnetic domain wall, which 
corresponds to a soliton solution of the nonlinear Landau- 
Lifshitz-Gilbert (LLG) that can preserve its shape during 
the propagation. Below we will see that all these results 
are modified when we introduce the magnetic field along 
the z-direction. 

Next, we apply a constant magnetic field h along the 
z-direction. By interpreting the dynamics of the spinon 
in terms of a real particle in a constant external field 
suggests the observation of Bloch oscillations. First we 
fix g = 0.3 and investigate the dependence of the spinon 
dynamics on h. The time evolution of the spinon under 
different magnetic fields h is shown in Fig. [3][a) h = O.I; 
(b) h = 0.2. For Bloch oscillations (BO) the frequency of 
the BO is only dependent on the strength of the external 
constant field, which can also be verified via our numer- 
ical result, as shown in Fig. |4l where we fix h = 0.2. We 
can find that for different g, the period of the BO is the 
same (T^ = 15.5), while the amplitude of BO A is deter- 
mined by g. If we fix 5 = 0.3, the time evolution of the 
spinon for different h is shown in Fig.lSJa). Following the 
interpretation in terms of a Bloch oscillation, we expect 



that the oscillation frequency to be proportional to the 
strength of the 'external field' h, which can be verified 
numerically (Fig. [5l[b)). 
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FIG. 6: The time evohition of the width A of the spinon for 
g = 0.1 and (a) h = 0.0; (b) h = 0.2. 

Next, we study the time evolution of width of the 
spinon. The width of the spinon A can be defined simi- 
larly to the definition of the width of a wavepacket. We 
introduce rii = 0.5 — \S^\ and A = ^ (ri)^ — (n^), , where 

(O) — . The result is shown in Fig. [51 Without 

the magnetic field along the z-direction, /i = 0, we ob- 
serve (Fig.|6](a)) that the width oscillates strongly at first 
whereas after some relaxation time, it grows linearly in 
time, which means the spinon would become wider and 
wider in the process of propagation of the spinon. The 
situation is different for h ^ 0, as shown in Fig 15] (b), 
where after some short-time behavior the width of the 
spinon oscillates around an average value instead of di- 
verging, which corresponds to a solitonic breather mode: 
not only is the position of the spinon confined to oscillate 
around the center of the lattice by the magnetic field, but 
its width oscillates as well around an average value. 
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FIG. 7: The time evolution of the position of the spinon for 
g = 0.1, h = 0.2 (solid line) and h ^ 0.2 + 0.05 cos(cjot) . 

The dynamics of the spinon is even more interesting 
when the system is driven by a time-dependent magnetic 
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field h(t). In our case, we choose a periodic driving mag- 
netic field: h{t) = ho + Scos{ujt). Classically, it is known 
that if the frequency of the driving force can match the 
intrinsic frequency of the system, the system can accumu- 
late vibrational energy and even a small periodic driving 
force can produce large amplitude oscillations, in another 
words, a resonance occurs. This phenomenon can also be 
observed in our quantum system, where the intrinsic fre- 
quency corresponds to the frequency of Bloch oscillations 
of the spinon. As shown in FiglTl the amplitude of the 
oscillation of spinon is drastically enhanced when the fre- 
quency of the driving potential is comparable to wq, the 
frequency of BO corresponding to h = h^. Highly non- 
trivial dynamics may emerge when the frequencies of the 
external driving force and the Bloch oscillation become 
incommensurate, which could induce chaotic dynamics of 
the spinon. 

Now we discuss several possible trapped ultracold 
atoms systems that could realize our Hamiltonian ([1]) ex- 
perimentally. Likely candidates are Rydberg atoms fll| . 
which have been already used in quantum simulations 
and manipulation and have attracted lot of attention in 
Albeit highly excited, the lifetime of 



12, He 



recent years 

Rydberg atoms can reach even hundreds of microseconds, 
which is much longer than the typical lifetime of the first 
excited state of alkali metal atoms. Two possible states of 
a Rydberg atom (excited state and ground state) on each 
site reduce the problem to that of a pseudo-spin system, 
and the strongly dipole-dipole interactions contribute to 
the interaction between the pseudo-spins (Jafcr|), and 
the external magnetic fields erf and erf can be realized by 
a laser with single atom Rabi frequency JIq and detuning 
A respectively. Notice that the dipole-dipole interaction 
drops off as r~^, meaning that the next-nearest neighbor 
interaction will still be appreciable, however it is much 
weaker than the NN neighbor interactions. Therefore, 
we do not expect that this changes the physics quali- 
tatively for a ferromagnetic interaction. An alternative 
way would be provided by ultracold dipolar atoms or 
molecules with two internal states in optical lattices in 
the hardcore limit, with an externally driven intra-species 
transition (g-field) and a Zeeman term[l3|. The recent 
experimental progress in observing ultracold atoms at the 
single atom level (or even single spin level) would 
be helpful to detect the dynamics of the single quasi- 
particles (spinous). The classic counterpart, a magnetic 
domain wall, has been realized experimentally as a super 
cold atom thermometer [3]. Considering the rapid devel- 
opment of the field, we expect Hamiltonian ([1} and the 
nontrivial spinon dynamics proposed in this paper should 
be accessible experimentally in various setups soon. 

Our results could also be illuminating for condensed 
matter physics, especially for the domain wall motion 
in magnetic nanowires. Classically, the magnetiza- 
tion dynamics of magnetic systems is governed by LLC 
eauation[l9j. where the spin is considered as a classic 



vector and its dynamics is reduced to a classic nonlinear 
equation. However, with advances in the miniaturization 
of magnetic structures one has to anticipate that the clas- 
sical and phenomenological LLC equation will eventually 
be inadequate, necessitating a full microscopic quantum 
description of the magnetization dynamics. Our result 
actually paves the way to study the motion of the do- 
main wall in the quantum situation. While dissipation 
would be very weak in a quantum optical implementa- 
tion of our model, it would be more relevant in a solid. 
For the classic domain wall, the dissipation is known to 
play a key role in determining its velocity [l^. Under the 
assumption of a memory-free bath, dissipation could be 
modeled in the framework of a Lindblad quantum mas- 
ter equation or a stochastic Schi^odinger equation [21|; 
the numerical method used here can be extended to sim- 
ulate this model, both in an matrix product operator 
(MPO)[22, 23 1 or quantum jump approach Depend- 
ing on the detailed nature of the physical realization (and 
hence the bath), much richer physics could emerge due to 
the interplay between interactions and dissipation mak- 
ing it an interesting topic for future studies. 
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